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I. INTRODUCTION 



The propulsion system of a ship is a complex, nonlinear system with 
several internal feedback loops. The performance of the ship itself is 
dependent, of course, on the capabilities and limitations of the pro- 
pulsion plant, which may limit overall ship performance so that the 
capabilities of the hull design are not realized. 

The first objective of this thesis is to model, in detail, a 
typical ship propulsion system, and to simulate it in the digital 
computer.'*' After validation of the computer model a simulation of the 
Mariner hull is to be coupled to the propulsion plant so that control 
system studies may be undertaken. As part of the validation procedure, 
behavior of the shaft speed was investigated with and without the 
governor . 

The second objective was to use the computer model as part of a ship 
feedback control system. For this purpose it was decided to explore the 
use of a feedback loop to maintain constant ship speed in turns. 



‘''IBM 360/67 DSL 360 LANGUAGE at W. R. Church Computer Center, Naval 
Postgraduate School. 
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II. PROPULSION PLANT DYNAMICS 



A. SHIP PROPULSION EQUATIONS 

The propeller angular speed is denoted by n. 

The basic equation describing the angular acceleration of a propul- 
sion shafting system is: 

"ZQ = 27rldn/dt 

Where: ZQ = Summation of torque - Qe-Qp 

I = Polar mass moment of inertia 
n = Propeller angular speed 
Qe = Shafting torque 
Qp = Propeller reaction torque 

then dn/dt = EQ 

2ttI 

For the typical gas turbine, the experimental data on the power 
turbine torque Qe versus the power turbine speed is shown in Fig. Al. 

The torque map of Fig. Al allows the engine torque to be determined if 
the engine speed and fuel flow rate Wf are know. However, this torque 
representation is correct only for steady state contitions. 

Propeller speed n(rps) is related to the engine speed N3(rpm) 
through the reduction gear ratio Kg 

N3 = 60Kg. n 

The block diagram for the ship propulsion system is given in Fig. A2. 

The power turbine speed can be controlled by a governor acting on 
the gas flow. The changes in fuel flow and speed that occur can be 
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assumed small enough to permit the use of constant coefficients in the 



dynamic equations. 

Using the time lag of the governor, t, the block diagram for the 
Governor is shown in Fig. A3, 



Where: 


n* = Rotational speed order 

n = Rotational speed 


t 


ner = n* -n = Speed error signal 

Kf = Fuel flow rate/RPM 

x = Time delay 

Wf = Fuel flow rate 

GOVERNOR = Kf/(l+xS) 



Then the propulsion plant is shown in Fig. A4. 
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Figure Al. Engine Torque versus Engine Speed 
and Fuel Flow Rate 



15 




Figure A2. Block Diagram of Ship Propulsion System 
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Figure A3. 




Block Diagram for the Governor. 
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GOVERNOR 




Figure A4. The Propulsion Plant with Governor 
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III. HULL EQUATION AND SIMULATION 



WAKE FRACTION : W 

The difference between the ship speed and the speed of advance of 
the propeller Vp is called Wake speed. 

The Wake fraction is defined as 
W - (V-Vp)/V 
thus Vp=V(l-W) 

THRUST COEFFICIENT : Ct 
TORQUE COEFFICIENT : Cq 

For the propeller in open water 



Ct 





( 1 ) 



Cq 



Q 

3 2 2 2 

pD (Vp +n D ) 



( 2 ) 



Where T = Thrust 
Q = Torque 

p = Water mass density 
Ct = Thrust coefficient 
Cq = Torque coefficient 
D = Propeller diameter 



From (1) and (2) 

2 2 2 2 

T - Ct.pD (Vp +n D ) 

3 2 2 2 

Q = Cq.pD (Vp +n D ) 
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The experimental data on the Wake fraction versus the speed of the 
ship is shown in Fig. Bl. 

The Thrust coefficient Ct and Torque coefficient Cq versus the 
Second modified advance coefficient are shown in Fig. B2 and Fig. B3. 
Where Vp = V(l-W) 

And a = (n.D) /Vp^+n^D^) 

Total ship Resistance, Rt, is a computer look-up table giving Rt 
versus Ship speed V (fig. B4) . 
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Figure Bl. Wake Fraction and Thrust Deduction 
versus Ship Speed 
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Figure B2. Thrusting Coefficient Ct versus Second 
Modified Advance Coefficient a 
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Figure B3. Torque Coefficient Cq versus Second Modified 
Advance Coefficient a. 
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Figure B4. 



Ship Resistance versus Ship Speed V. 
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IV. COMBINED PROPULSION PLANT AND HULL EQUATION 



A. POWER TURBINE WITHOUT GOVERNOR : 

Nonlinear functions such as Wake fraction, Thrust deduction fraction 
Thrust coefficient, torque coefficient, ship resistance and engine torque 
were stored in look-up tables during the computer solution. 

Using the DSL/360 language these equations were solved by means of 
an integration and updating from table look-up every 0.1 sec until the 
steady state was reached. 

Fig. B5. shows propeller shaft speed (rpm) versus time 
Fig. B6. shows ship speed V versus time 
Fig. B7. shows propeller action torque versus time 
Fig. B8. shows shaft torque Qe versus time 
Fig. B9. shows Cq versus time 
Fig. BIO. shows Ct versus time 
For a step change in. fuel flow rate 
Wf = 7000. + 3000. step (0.0) 

The shaft torque decreases and goes to steady state. The propeller 
action torque increases, overshoots at 10 sec. then goes to steady state. 
The propeller shaft speed increases, goes to steady state in about 80 sec 
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Figure B5. Shaft Propeller Speed N4 versus Time. 
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Figure B6. Ship Speed versus Time. 
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Figure B7. Shaft Torque versus Time. 
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Figure B8. Propeller Reaction Torque Qp versus Time. 
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Figure B9. Thrust Coefficient Ct versus Time. 
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Figure BIO. Torque Coefficient Cq versus Time 
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B. POWER TURBINE GOVERNING: 



Turbine speeds with governor are shown in Fig. A4 and compared with 
the responses that occur without a governor. 



For small step changes in 


power demand at 7000 Lb/hr fuel flow rate. 


the governor, which has a rate 


of 40 Lb/hr. RPM and time lag 1/^ =1/5=0. 2 


sec., causes the power turbine 


response to be reduced from 80 sec. to 


20 sec. with an overshoot of 3. 


,92% of the final change in speed. 


RPM command 





n* = 120. +4. STEP (0 . ) 



Now changing 




A = 45. 


. RPM (Fig. Cl) 


o 

(1 

<] 


, RPM (Fig. C2) 


t> 

ii 

u> 

o 


. RPM (Fig. C3) 


A = 20. 


, RPM (Fig. C4) 


A = 10. 


, RPM (Fig. C5) 



Percent Overshoot as a function of rudder angle is shown on Table 1. 
The block diagram for the combined propulsion plant and hull 
equations were shown in Appendix B. 
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TABLE 1. Percent Overshoot with Governing. 
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Figure Cl. Shaft-Speed versus Time, with 45 RPM Increase 
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Figure C2. Shaft Speed versus Time, with 40 RPM Increase 



CpOl 

f ' » 1] 

u i 




j 



i 



j 

I 

1 



I 






I 




o 

cr 



o 

o 



o 

o 

rO 



CD 

Cn 



CD 

CD 

ID 



CD 

CD 



CD 

CD 

<N» 



36 



Figure C3. Shaft-Speed versus Time, with 30 RPM Increase 
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Figure C4. Shaft-Speed versus Time, with 20 RPM Increase 
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Figure C5. Shaft-Speed versus Time, with 10 RPM Increase 



V. EQUATIONS OF MOTION OF THE HULL 



A. NEWTON’S LAW OF MOTION: 

Newton's Law of motion for a rigid body can be written as two 
equations: 

FORCE EQUATION 

^ v 

F = Force on body = d (Momentum) /dt 
= d(mU^) 

Where U„ Velocity of body 

Lr 

MOMENT EQUATION 

M = Moment acting on a body 

d (Angular momentum) 
dt 

= d(m) 

dt 

Where I = Moment of inertia 
.52. = Angular velocity 



B. SIX DEGREES OF FREEDOM: 
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Redrawing: 




Z 

Y 

X 

4 > 

e 



Vertical axis, positive downward 
Tranverse axis, positive staboard 
Longitudinal axis, positive forward 
Roll angle 



Positive rotations are indicated 



Pitch angle 
Yaw angle 



in the sketch above 



t = = mU + X(mU) 

at 

M = - d -y H ^ = H + SI x H 
dt 



m: Mass of ship 
U = ui+vj4vk (3) 

u: Rate of surging 
v: Rate of swaying 
w: Rate of heaving 
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R=pi+qj+rk 



(4) 



R: Angular velocity 

p = = Rate of roll 

q = 0 = Rate of yaw 
r = 4^ = Rate of pitch 



F = Xi + j + Zk 



(5) 



X: Hydrodynamic force along X axis 
Y : Hydrodynamic force along Y axis 
Z: Hydrodynamic force along Z axis 



M = Moment vector acting on body 
= Ki + Mj + Nk (6) 



K: Rolling moment about X axis 

M: Pitching moment about Y axis 

N: Yawing moment about Z axis. 



Substituting equations (3) , (4) , (5) , and (6) into equation (2) 

gives 

r\ a 

X = m(u + qw - rv - X G (q + r ) + Y G (pq - r) + Z G (pr + q)) 

^ Y = m(v + ru - pw - Y G (p 2 + r 2 ) + Z G (qr - p) + X G (qp + r)) 

Z = m(w + pv - qu - Z G (p 2 + q 2 ) + X G (rp - q) + Y G (rq + p))' 

K = I x p + (I z - I y ) qr - m(Y G (w + pv -qu) - + ru - pw)) 

M = I y q + (I x - I z ) rp - m(Z G (u + qw - rv) - X Q (w + pv - qu) ) 

N = I z r + (1^ - I x ) pq - m(X G (v + ru - pw) - Y G (u + qw - rv) ) 



For the surface ship in calm water, roll, pitching and heaving are 
taken as zero. That means p=q=w=0 and most ships have Y^= 0. 
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Then equation (7) becomes: 



X = m(u-rv-X G r 2 ) 

Y = m(v+ru+X r ) (8) 

C? 

N = I 2 r-m(X G (v+ur)) 



C. LINEARIZATION OF THE EQUATIONS OF MOTION: 

1. X Equation 

The nonlinear TAYLOR expansion for the X equation is 



X r X°+ + 

-f — (Aty t Xw ^ + Yii ^ ^ 

+ 2 Yu, v u . • ♦ • + 2. Xf.s 



(5) 



«2 

+ s x uuv t i W‘** + •••+ ** 

+ + j 



Note: 



1. The TAYLOR expansion of the equation above includes terms up 

to third order, since terms higher than third order are not significant. 

2. Since X is an EVEN function of v, of , of 

X(v) = ^ + vV &6 V + - - * 

X(r) = ^ t ^ 

x(i) = Cj S l + C H s\ C& /+ - • • 



Then the odd powers of v,r,<f are zero i.e. all X^VjX^r ,X<f <f are equal 
zero. 



42 



Thus equation (5) becomes: 



+x u ,W'> + Xu **(6) 

Some of these terms are equal to zero (see Table 2) 
then X=X Q =0 



2. Y Equation 



The nonlinear TAYLOR expansion for the Y equation is: 



y= y" + [y«au +xv+y*'* +>}<*] 

+J. f* y (4u) + X'i* ^ ^ 'i’vfl' 4 4 ' ' * 

/MR 

+ 2)^*0 + 




y uw 4 y v ^ v3 + 

+ 3 y«nv (a “' )lv + S + ' --1- } hs if 

+ 6 y J ( 1 ) 



NOTE: 

1. The terms higher than third order are not significant. 

2. Y is an odd function of v,r >S 

Y (v) = dj $ 4 ^ -f ^ 5 * - * • 

Y( r ) = e, n + £>> rt+ % * * * 

y ( ) = 4 £ £**+ • * * 
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Then equation (7) becomes: 



Y- y° + 1u + /✓ V 4-y \n + y^ j + 

^ y^uu ^ /m v + ^ 

+ i >U H V + i +i IvAptJ’f 

+ 4- %u vf+s a*« « A w ^ 



Where the indicated terms are equal to zero (see Table 3) 

i--% f >^+y** +/,>* + £v w A 

( 8 ) 

3. N Equation 



The nonlinear TAYLOR expansion for the N equation is: 



Nr: A/ A/14 + Vy ^ + N T ft -f- f\l$ b + 



1 4- U/„„ («r+ u w * t u w * + Nu i'4 

o( L - 

+ <3 + •■ ■■ +c2 NtjiS j 

+•■“ fruuu AU * fAW^ + H SgS 6 h 
4 3 NuufW** 4 3 diiutM**-* . • • 

•f g Mji/Jfc, ^W. 4 - 



4 bNiffiS 1 

( 9 ) 



NOTE: 



1. The terms higher than third order are not significant. 

2. N is an odd function of v,r , S 

N(v) = 4 %b 1/5 4 %^ Sjt ‘ ' ‘ 

N(r) = H- + ^4 &5 ^ '"*■ 

N(«) = K> & 4 + M*+ ' ' 
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Then equation (9) becomes: 

hlz N u Au + N v v + N z n, 4 - + 

Z Huvu t £ Nnv ^ + i-/W *? + £ %| ^ 

+ J-^UUV ^ V+^J-Nhhs du ' 2 t ■ • ' - f N^S'/nS 
The indicated terms are equal to zero (see Table 4) 

Then 

N z AfirVf Hjt +Nii + A- Nvw *N- g- Mrtf't 4 - + J 

( 10 ) 

Combining the third order TAYLOR expansion for X,Y, and N with the 
dynamic response terms of the X,y, and N equation. 

X equation X = 0 
Y equation 

(m-Yv)v+(mX G -Yp r = Y(u,v,r,<0 

N equation (1) 

(mX G -N^)v+(I z -Npr = N(u,v,r,<T) 
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Using CRAMER’S rule to solve these equations 



Y(u,v,r,<f) “Xq-Y,- 

N(u,v,r ,£) I z -N r 



m-Y- 

V 



mx G -Y 



r 



mX -N- 
G v 



I--N- 
Z r 



m-Y^ Y(u,v,r, ) 

mX G -N^ N(u,v,r, ) 



r = 



m-Y- 

v 



mX G - 



Y- 

r 



mX -N- 
G v 



I z -N 



r 
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TABLE 2. ASSESSMENT OF THE COEFFICIENTS IN 
THE X EQUATION 



Variable 


Coefficient 


Series 60/Model 5.1.1 


ix 


m-Xu 


0.1750 




Au 


Xu 






flu 2 


1/2 Xuu 






flu 5 


1/6 Xuuu 


' 




V 2 


1/2 Xw 






r 2 


1/2 Xrr + mX„ 








1/2 Xj* 






V V 2 A U 


1/2 Xvvu 






Au . 


1/2 Xrru 






,2 

S A U 


1/2 X#u 






v /r 


Xvr+m 






v<r , 


Xvj 






rj 


X r<T 






vr^u 


Xvru 






V^AU 


Xv u 






r/au 


Xr u 






— 


Xo 







Note: No entry in these columns means the coefficient was ignored. 



*From PRINCIPLES OF NAVAL ARCHITECTURE, SNAME, 1968, page 548. 
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TABLE 3. ASSEMSSEMNT OF THE COEFFICIENTS IN 
THE Y EQUATION 



Variable 


coefficient 


Series 60/5.1. 


• 

V 


m-Yv 


0.509 


Y9 


• 

r 


mX -Yr 




Y8 


V 


Yv 


-0.260 


Y1 


v 3 


1/6 Ywv 


-2.150 


Y2 


vr^ 


1/2 Yvrr 


-1.180 


Y3 


vS 2 


1/2 Y v« 






VAU 


Yvu 






. 2 

VAU 


l/2Yvuu 






r 


(Yr-m) 


-0.0781 


YY4 


r 3 


1/6 Yrrr 


-0.0461 


Y5 


TV 2 


1/2 Yrvv 


-0.0994 


Y6 


rS 2 


1/2 Yr SS 






rAu 


Yru 






i*a u x 


Yruw 






6 


1/6 % S 




Y7 


6 v 2 


1/2 Y«frr 








Y«Tu 






S AXL 2 


1/2 YJuu 






v r<T 


YvrJ 






— 


Yo 


0.00016 


Yo 


AM, 


YY’ 






(Auf 









*No entry in these columns means the coefficient was ignored. 
*From PRINCIPLES OF NAVAL ARCHITECTURE, SNAME, 1968, page 548. 
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TABLE 4. ASSESSMENT OF THE COEFFICIENT IN 
THE N EQUATION. 



Variable 


coefficient 


Series 60/Model 5.1.1 


V 


mX r -Nv 




N9 


• 

r 


i*-»* 


0.0202 


N8 


V 


Nv 


-0.0750 


N1 




1/6 Nwv 


-0.5850 


N2 


V r 2 
v/ 2 
Vdu 


1/2 Nvrr 
1/2 Nv# 
Nvu 


-0.5060 


N5 






VAU^ 


1/2 Nvuu 






r 


(Nr-mXg) 
1/6 Nrrr 


-0.0569 

-0.1010 


NN4 

N5 


rv 2 

IS 2 

1'6U 

2 

tau 


1/2 Hrw 
1/2 Nr<tf 
Nru 

1/2 Nruu 


-1.4200 


N6 










S 


V$ 


-0.0240 


N7 


o 


1/6 








1/2 N£vv 
1/2 N<5rr 
N£u 

1/2 N<fuu 
Nvr<$ 






<M 

in 






Ja u 






2 






vr$ 






— 


No 


-0.0005 


No 


4U 


NS 






c ^«> 1 


nSh 







*No entry in these columns means the coefficient was ignored. 
*From PRINCIPLES OF NAVAL ARCHITECTURE, SNAME, 1968, page 549. 
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D. RELATIONS BETWEEN FIXED AXES AND MOVING AXES: 



So. far the reference axes X,Y, and Z are fixed to the moving ship. 

In order to determine the path of the ship, it is convenient to take the 
fixed axes to the earth X^,Y^. 

In figure Dl., the ship is at the position A. 

The YAW ANGLE ^ is the angle between the U axis and the tangent line 
with the path at point A. 

Then the following relationships between the fixed and the moving 
axes can be written: 

X 0 (t) = U(t) Cos t(t)-V(t)Sin^ (t) 

Y q ( t) = U(t) Sin l / / (t)+V(t)Cos '/'(t) 




Figure Dl. Relations Between Fixed Axes in the 
Earth and the Ship Axes. 
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E. RUDDER DEFLECTION: 



Let S - Rudder angle 

DRATE = Rudder deflection rate. 

TDMAX = Time at maximum rudder deflection 



Then the rudder deflection in turning maneuver can be shown in 



Figure D2: 




Let Rudder deflection = 8 

DRATE = Rudder deflection rate = D /TDMAX 

max 

Then the function of the rudder deflection can be written: 



S = DRATE Tramp (t. )-RAMP(t-t )| -r^-r 
I lag max J 57.3 
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F. THE TURNING PATH OF A SHIP: 



There are 4 phases of a turn: 




Figure D3. The Turning Path of a Ship. 

Approach phase V = 0 , r = 0, v = 0, r = 0 

• • 

First phase of turn v/0, r ^ 0, v = 0, r = 0 
Second phase of turn v ^ 0, r ^ 0, v ^ 0, r/0 

Third phase of turn v = 0, r = 0, v ^ 0, r^O 

Using dimentional equation (1) and introducing the rudder deflection 
30°. The characteristics of transient phases of a turn are shown in 
Figures D6, D7, D8, and D9 . 

• • 

A-B Second phase v / 0, v ^ 0, r^0, r^0 

B-C Third phase v ^ 0, v = 0, r^0, r^0 

Now using different rudder deflections =35°, 30° the turning path 
of a ship is shown in Figures D4 and D5. When using the bigger deflection 
the steady turning radius is smaller. 
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Figure D4. The Turning Path of a Ship with 35 Degree Rudder Angle. 
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Figure D5. The Turning Path of a Ship with 30 Degree Rudder Angle. 
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Figure D6 • Angular Acceleration about Z Axis. 
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TURN 

R 




Figure D7» Angular Velocity about Z Axis. 
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Figure D8. Angular Acceleration about Y Axis. 
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VI. SPEED CONTROL IN TURNS 



A. DYNAMIC BEHAVIOR OF THE PROPULSION PLANT: 

The propulsion plant and turning path were observed. Now combine 
both of them and add feed-back of ship speed with Gain K1 in order to 
keep the ship speed constant in turns. 



The block diagram for this part can be drawn: 




Figure El. Feedback Control Velocity. 



Where : 



V*: Reference ship speed. 

V : Ship speed. 

n*: Shaft angular speed command. 

U : Velocity about X axis. 

VA: Velocity about Y axis' 

Kl: Ship speed gain 

n : Shaft angular speed 

6 : Rudder angle 



The block diagram of this system was shown in Appendix B. 
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B. SIMULATION METHODS 



1. Determine all initial conditions in the steady state 

RUN 1. 

Let propeller shaft speed command n* be constant 
n* = 130 RPM 

And no feedback of ship speed 
K1 = 0 

Then determine the steady state speeds (i.e., determine 
the component of initial condition ship speed on X axis 
and on Y axis (UIC, VIC). 

RUN 2. 

Now the propeller shaft speed command n*-130 rpm making 
K1 = 8,16,24,32 respectively. Then determine the steady 
state in each case (i.e., find all IC’s in the steady 
state) . This guarantees that the loop is in steady 
state, since the initial conditions on other variables 
might still cause a transient. 

2. Put in a step change in the ship speed with 2 kts. increase and 
delay time 10 sec. Then V* becomes: 

V* = 20. + 2.0 step (10.) 

C. SPEED CONTROL STUDIES 

Using the inital conditions of RUN 1 and 2, a rudder command was 
applied at t = 0, causing a speed decrease due to added drag. At 
t = 10 sec, the commanded speed, V* was increased by 2 knots, stepwise. 
This caused a speed increase. The effects of changing the loop gain, 

Kl, is shown in Figure E2. 
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D. DISCUSS THE CURVES 

1. Ship resistance Rt, Wake fraction W, Thrust fraction t 

The ship speed change was small, so that: 

a. The ship resistance Rt = K.V (see Fig, B4) K = constant 
then the form of Rt was the same as V (Fig, E3 curve 1) 

b. The Wake fraction was constant (Fig. Bl) 

then W = (Fig. E3 curve 2) 

c. In this interval the slope of the thrust fraction was 
negative. Then when V decreases to increases and as V 
increases t decreases (Fig. E3 curve 3) 

2. o curve, torque coefficient Cq, thrusting coefficient 

VP = (1-W) V 

= (1 -.0*5) V 

= ^55M 

nD 

and a = — ~ 

Vp -n D 

The a curve is shown in Fig. E4 curve 1 

a changed from.^to .320 then the slope of Ct and Cq curve was constant. 
Then the Ct and Cq curve have the same form as the <T curve. 

3. Propeller action torque Qp, shaft torque Qe and propeller angular 
speed n 

3 2 2 2 

Qp = Cq.p.D (VP +n D ) 

The Qp curve was shown in Fig. E5 curve 2. 

The Wf curve and n curve are shown in Figures E6 and E7. 
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Figure E2. Ship Speed V versus Time, with K1 = 8,16,24,32 
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Figure E3. Ship resistance, Wake Fraction and Thrust Fraction versus Time. 
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Figure E4. 



Second Modified Advance Coefficient , Torque Coefficient 
Cq, Thrust Coefficient Ct versus Time. 



65 




Figure E5. Propeller Action Torque Qp and Shaft Torque versus Time. 
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Figure E6. Fuel Flow Rate Wf versus Time. 
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Figure E7 . Propeller Angular Speed versus Time. 
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Figure E8. NER versus Time 
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Figure E9. U, VA, R versus Time 



VII. CONCLUSION 



This thesis has presented the combined model of a propulsion plant 
and hull. DSL/360 LANGUAGE was found to be powerful for simulating ship 
motion. 

In turbine powered ships, it is necessary to predict the propeller 
response to changes in power demand. 

This problem has been investigated but the studies have been con- 
centrated on the effect of dynamics of the ship. This work is concerned 
both with the dynamic of the propulsion plant and ship dynamics of the 
ship. 

It has been shown that the complex, nonlinear dynamics of propulsion 
plant and hull can be adequately represented by an all digital computer 
model. The model permits studies of the internal dynamics of the pro- 
pulsion plant, as well as the external dynamics of the hull. 

It has been shown that the model can be incorporated in a feedback 
control system and its effects included in studies of control system 
performance. 
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VIII . RECOMMENDATIONS 



The combination of propulsion plant and hull in a single computer 
model provides a realistic tool for many studies. The following are 
topics suitable for future investigation. 

a. Course keeping, station keeping and replenishment at sea. 

b. Effect of propulsion plant dynamics on maneuvering in sea 
states (or regular waves) . 

c. Effect of sea state conditions on the operation of the 
propulsion plant. 

d. Discuss variation in steady state speed with increase command 
as shown in Fig. Cl - C5. Also discuss effect of governor 
loop on steady state speed. Fig. E2. 



72 



APPENDIX A 



Table Look-up for Ship Resistance versus Speed : 




DIMENSION VT(13) ,RTT(13) 

DATA VT/0 .0,2.5,5.0,7.5,10. , 12 .5 , 15 . , 17. 5 , 20. , 22 .5 , 25 .5, 127 . 5, 30. / 
DATA RTT/2*0. 0,7000. , 13000 . , 25000. , 39000. ,57000. ,80000. ,1103000., 
136000. ,173000. ,225000. ,280000./ 

IF(V.GT. 30) GO TO 1 
DELV 2.5 

N IFIX(V/DELV) 1 

SLOPR(RTT (N l)-RTT(N)) /(VT(N-l)-VT(N)) 

TABR SLOPR* (V-VT (N) ) -RTT (N) 

RETURN 

1 TABR 280000. 

RETURN 

END 
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CALL EMDRW ( NPLOT ) 
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